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We present a systematic study on the influence of epitaxial strain and hole concentration on the 
magnetic anisotropy in (Ga,Mn)As at 4.2 K. The strain was gradually varied over a wide range 
from tensile to compressive by growing a series of (Ga,Mn)As layers with 5% Mn on relaxed graded 
(In,Ga) As/GaAs templates with different In concentration. The hole density, the Curie temperature, 
and the relaxed lattice constant of the as-grown and annealed (Ga,Mn)As layers turned out to be 
essentially unaffected by the strain. Angle-dependent magnetotransport measurements performed 
at different magnetic field strengths were used to probe the magnetic anisotropy. The measurements 
reveal a pronounced linear dependence of the uniaxial out-of-plane anisotropy on both strain and 
hole density. Whereas the uniaxial and cubic in-plane anisotropics are nearly constant, the cubic 
out-of-plane anisotropy changes sign when the magnetic easy axis flips from in-plane to out-of- 
plane. The experimental results for the magnetic anisotropy are quantitatively compared with 
calculations of the free energy based on a mean-field Zener model. An almost perfect agreement 
between experiment and theory is found for the uniaxial out-of-plane and cubic in-plane anisotropy 
parameters of the as-grown samples. In addition, magnetostriction constants are derived from the 
anisotropy data. 

PACS numbers: 75.50.Pp, 75.30.Gw, 75.47.-m, 61.05.-a 
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I. INTRODUCTION 

Spin-related phenomena in semiconductors, such as 
spin polarization, magnetic anisotropy (MA), and 
anisotropic magnetoresistance (AMR), open up new con- 
cepts for information processing and storage beyond con- 
ventional electronics jiiSi Being compatible with the stan- 
dard semiconductor GaAs, the dilute magnetic semicon- 
ductor (Ga,Mn)As has proven to be an ideal playground 
for studying future spintronic applications. "^'^ In particu- 
lar, the pronounced MA and AMR, largely arising from 
the spin-orbit coupling in the valence bandj^i^ potentially 
apply in novel non- volatile memories and magnetic-field- 
sensitive devices. Ferromagnetism is implemented in 
(Ga,Mn)As by incorporating high concentrations (^1%) 
of magnetic Mn^+ ions into the Ga sublattice. The fer- 
romagnetic coupling between the S=5/2 Mn spins is me- 
diated by itinerant holes provided by the Mn acceptor 
itself. Curie temperatures Tq up to 185 K, i.e. well 
above the liquid-N2 temperature, have been reported^i^ 
and there is no evidence for a fundamental limit to higher 
values.^ 

The magnetic properties of (Ga,Mn)As are strongly 
temperature dependent and can be manipulated to a 
great extent by doping, material composition, and strain. 
In (Ga,Mn)As grown on GaAs substrates, however, hole 
density p, Mn concentration x, and strain e are intimately 



linked to each other and cannot be tuned independently 
by simply varying the growth parameters, p and e sen- 
sitively depend on the concentration and distribution of 
the Mn atoms which are incorporated both on Ga lattice 
sites (MuQa) and, to a lower extent, on interstitial sites 
(Mni), where they act as compensating double donors. 
Post-growth treatment techniques such as annealing or 
hydrogenation are frequently used to increase or decrease 
the hole concentration due to outdiffusion and/or rear- 
rangement of MnjiSiiiii^ii^ii^ Or duc to the formation 
of electrically inactive (Mn,H) complexes ; ^^d^d"^ respec- 
tively. In both cases, however, the treatment concur- 
rently leads to a decrease or increase of the lattice pa- 
rameter, respectively, and thus to a change of the strain. 

The epitaxial strain in the (Ga,Mn)As layers, aris- 
ing from the lattice mismatch between layer and sub- 
strate, can be adjusted by tailoring the lattice parame- 
ter of the substrate. While (Ga,Mn)As grown on GaAs 
is under compressive strain, tensily strained (Ga,Mn)As 
can be obtained by using appropriate (In,Ga)As/GaAs 
templates i^ d^d^i^o Experimental studies addressing this 
issue, however, have so far been restricted to merely a 
limited number of representative samples. 

In this work, the influence of epitaxial strain and hole 
concentration on the MA at 4.2 K is analyzed in a sys- 
tematic way by investigating a set of (Ga,Mn)As layers 
grown on relaxed (In, Ga) As/GaAs templates with differ- 
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ent In concentration. Keeping the Mn content at ~5% 
and changing the maximum In content in the (In,Ga)As 
buffer layers from 0% to 12%, the vertical strain ezz in 
the as-grown (Ga,Mn)As layers could be gradually varied 
over a wide range from £^2=0.22% in the most compres- 
sively strained sample to e^z^— 0.38% in the most tensily 
strained sample without substantially changing p. Post- 
growth annealing leads to an increase in p, yielding a 
second series of samples with nearly the same range of 
Ezz but higher hole concentrations. The strain depen- 
dence of the anisotropy parameters for the as-grown and 
the annealed samples was determined by means of angle- 
dependent magnetotransport measurements i^i^ Part of 
the experimental data has already been published in con- 
ference proceedings.— Here, we combine the earlier with 
the present extensive experimental findings advancing a 
quantitative comparison of the intrinsic anisotropy pa- 
rameters with model calculations for the MA, performed 
within the mean-field Zener model introduced by Dietl et 
ali^ Note that several samples analyzed in Ref. [2^ have 
been substituted by new samples grown under optimized 
conditions and that the sample series has been expanded 
by one specimen with e22=-0.38%. 



II. EXPERIMENTAL DETAILS 

A set of differently strained (Ga,Mn)As layers with 
constant Mn concentration of and thickness of 

'^ISO nm was grown by low-temperature molecular-beam 
epitaxy (LT-MBE) on (In,Ga)As/GaAs templates with 
different In content in a RIBER 32 MBE machine. 
Indium-mounted semi-insulating VGF GaAs(OOl) wafers 
were used as substrates. After thermal deoxidation, a 
30-nm-thick GaAs buffer layer was deposited at a sub- 
strate temperature of TgRiSSO °C. Then the growth was 
interrupted, Tg was lowered to ~430 °C, and a graded 
(In,Ga)As buffer with a total thickness of up to ~5 /im 
was grown. Starting with Ino.02Gao.98As, the temper- 
ature of the In cell was first continuously raised to in- 
crease the In content up to a value of <12% and was then 
kept constant until the required thickness of the buffer 
layer was reached. The growth was again interrupted, 
Ts was lowered to --250 °C, and the (Ga,Mn)As layer 
was grown in As4 mode at a growth rate of ~200 nm/h. 
The growth was monitored by reflection high-energy elec- 
tron diffraction showing no indication of a second-phase 
formation. The use of a graded (In,Ga)As buffer— min- 
imizes the deterioration of the (Ga,Mn)As layer caused 
by threading dislocations in the relaxed (In,Ga)As/GaAs 
template. The resulting (Ga,Mn)As layers exhibit nearly 
the same quality as conventional samples directly grown 
on GaAs.'^'^ After the growth, the samples were cleaved 
into several pieces and some of the pieces were annealed 
in air for 1 h at 250 °C. The structural properties of 
the (Ga,Mn)As layers were analyzed by means of high- 
resolution x-ray diffraction (HRXRD) measurements per- 
formed with a Siemens D5000IIR x-ray diffractometer us- 



ing the Cu-Kqj radiation at 0.154 nm. Hall bars with 
current directions along the [100] and [110] crystallo- 
graphic axes were prepared from the samples by stan- 
dard photolithography and wet chemical etching. The 
width of the Hall bars is 0.3 mm and the longitudinal 
voltage probes are separated by 1 mm. The hole densities 
were determined by high-field magnetotransport mea- 
surements (up to 14.5 T) at 4.2 K using an Oxford SMD 
10/15/9 VS hquid-hehum cryostat with superconducting 
coils. The Curie temperatures were estimated from the 
peak positions of the temperature-dependent sheet re- 
sistivities at 10 mTi^i2^iS£ The MA of the samples was 
probed by means of angle-dependent magnetotransport 
measurements at 4.2 K using a liquid-He bath cryostat 
equipped with a rotatable sample holder and a standard 
LakeShore electromagnet system with a maximum field 
strength of 0.68 T. To determine the saturation magne- 
tization, we employed a Quantum Design MPMS-XL-7 
superconducting quantum interference device (SQUID) 
magnetometer using the Reciprocating Sample Option 
(RSO). The measured SQUID curves were corrected for 
the diamagnetic contribution of the substrate. 

As discussed in detail in Ref. 113, our samples ex- 
hibit spin wave resonances which are most pronounced 
for the external magnetic field oriented perpendicular 
to the sample plane. These spin wave excitations have 
been traced back to an inhomogeneous free-energy den- 
sity profile, or more precisely to a linear variation of the 
MA parameters along the growth direction, presumably 
arising from a vertical gradient in the hole density.— 1^ 
Therefore, all physical parameters derived via magneto- 
transport in this study have to be considered as ejfec- 
tive parameters representing the averaged electronic and 
magnetic properties of the layers. 



III. THEORETICAL CONSIDERATIONS 

In the present context, MA represents the dependence 
of the free-energy density F on the orientation m of the 
magnetization M = Mm^ In the absence of an ex- 
ternal magnetic field, m is determined by the minimum 
of the free energy. Since for reasons of crystal symme- 
try F usually exhibits several equivalent minima, more 
than one stable orientation of M exists. This symmetry- 
induced degeneracy of F can be lifted by the application 
of an external magnetic field H. 

The theoretical considerations on the AMR and the 
MA in this paper are based on a single-domain model 
with a uniform magnetization M. While the direction 
of M is controlled by the interplay of F and H, its 
magnitude M is assumend to be constant under the 
given experimental conditions. For sufficiently high field 
strengths H, this assumption can be considered as a 
good approximation. Hence, the normalized quantity 
Fm=F/M is considered instead of F, allowing for a con- 
cise description of the MA. 
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A. Phenomenological description of the MA 

There are several contributions to Fm which we refer 
to as intrinsic (magnetocrystalline) or extrinsic: 

Pm — FM,int + FM,oxt ■ (1) 

The intrinsic part FM,int=FM,c+FM,s originates from the 
holes in the valence band {Fm,c) (see Sec. lIIIB]) and from 
the localized Mn spins (Fm.s)'^ Whereas Fm.c is strongly 
anisotropic with respect to the magnetization orienta- 
tion, reflecting the anisotropy of the valence band, the 
localized-spin contribution Fm,s=Jo dM'noH{M')/M 
is isotropic and therefore irrelevant for the following con- 
siderations. In a phenomenological description, FM,int 
can be expressed in terms of a series expansion in ascend- 
ing powers of the direction cosines TOj., rriy, and to^ of 
the magnetization with respect to the cubic axes [100], 
[010], and [001], respectively. Considering terms up to 
the fourth order in m, i^M.int for cubic systems with 
tetragonal distortion along the [001] growth direction is 
given by^^ 

FuM-m) = Bo+B2i_ml+B^{ml+ml)+Bi^m'l . (2) 

In the case of a perfect cubic crystal, symmetry requires 
i32_L=0 and B^—B^^^. 

The extrinsic part i^M.oxt comprises the demagneti- 
zation energy due to shape anisotropy and a uniaxial 
in-plane anisotropy along [1 10]. The origin of the lat- 
ter anisotropy is controversially discussed. It is traced 
back either to highly hole-concentrated (Ga,Mn)As clus- 
ters formed during the growth)^ to the anisotropy of 
the reconstructed initial GaAs (001) substrate surfaccj^i 
or to a trigonal-like distortion which may result from a 
nonisotropic Mn distribution, caused, for instance, by the 
presence of surface dimers oriented along [110] during the 
epitaxy 42 Approximating the (Ga,Mn)As layer by an in- 
finite plane, we write the total extrinsic contribution as 

FM,ext{rn) = Bdml + BiiJ^{m.^ - niy f, (3) 

where Bd=ii[)M /2. 

In the presence of an external magnetic field H—Hh, 
the normalized Zeeman energy —^oHm has to be added 
to the total free-energy density. This corresponds to a 
transition from Fm to the normalized free-enthalpy den- 
sity 

-Booi 

GA/(m) = Bo + {B2± + Bd)ml 

+B4\\{ml + wiy) + B^j^ml 

+BiiQ^{mx - niy)^ - fioHhm. (4) 

The anisotropy parameters i?2± and B^ are both related 
to rnj, and are therefore combined into a single parameter 
-Bqoi. Given an arbitrary magnitude and orientation of 



H, the direction of m is determined by the minimum of 
Gm- 

All anisotropy parameters introduced above are in 
SI units. Expressed by the aniso trop y fields in 
cgs units as defined, e.g., in Ref. [SJ, they read 

as Biio=-/ioH2||/2, B2±=-fiQH2±/2, B4||=-^o-f^4||/4, 
and B4±=—iJ.oH4±/4:. Note also that the magnetic 
anisotropy field used in Ref. ^23 and the anisotropy pa- 
rameters used here are related via /zo-ffaniso — '^i^cS^ + 
K^,)/M,,t^-~2Booi~4B4±. 



B. Microscopic theory 

For a microscopic description of the intrinsic part 
Fm ;mt, we adopt the mean- field Zener model of Dietl et 
al. introduced in Ref . i5|. The objective of the microscopic 
calculations discussed below is first, to justify the approx- 
imation in Eq. made by considering only terms up 
to the fourth order, and second, to compare the experi- 
mentally found dependence of the intrinsic anisotropy pa- 
rameters B2±, -B411, and B4± on e^z andp (see Sec. lIVC]) 
with that predicted by the mean-field Zener model. 

According to the k ■ p effective Hamiltonian theory 
presented in Ref. |^, the Hamiltonian of the system is 
given by 

H = Hkl + Ti-e + Hpd- (5) 

Here, Hkl represents the 6x6 Kohn-Luttinger k p 
Hamiltonian for the valence band and He=X]i j D'^^^^Sij 
accounts for the strain Sij in the (Ga,Mn)As layer via the 
deformation potential operator Tipd=—No(3Ss de- 

scribes the p-d hybridization of the p-like holes and the 
localized Mn d-shell electrons, which results in an in- 
teraction between the hole spin s and the Mn spin S 
carrying a magnetic moment SgfiB- Here, g=2 is the 
Lande factor, /ie the Bohr magneton, and f3 and Nq de- 
note the p-d exchange integral and the concentration of 
cation sites, respectively. In terms of the virtual crystal 
and mean-field approximation, the exchange interaction 
can be written as Ti.pd~Msf3/gfiB- Explicit expressions 
for the individual contributions in Eq. ([5|) can be found in 
Ref. i5|. As an approximation, the values of the Luttinger 
parameters 7^ {i— 1,2,3), the spin-orbit splitting Aq, and 
the valence band shear deformation potential b are chosen 
as those of GaAs. Explicit values are 71=6. 85, 72=2.1, 
73=2.9, Ao=0.34eV and 6=-1.7eV, respectively,^ The 
quantity parameterizing the exchange splitting of the va- 
lence subbands is given by 

Bg = -z , (6) 

with the Fermi liquid parameter A-p. In contrast to 
Ref. i, we restrict our calculations to zero temperature 
(T=0) and zero magnetic field (H=0). In this approxi- 
mation, the Fermi distribution is represented by a step 



4 



function and the Zeeman as well as the Landau split- 
ting can be neglected. These simplifications are justified, 
as our measurements were carried out at T—4:.2 K and 
fj-oH < 0.7 T, where the Zeeman and Landau splittings 
of the valence band are expected to be much smaller than 
the splitting caused by the p-d exchange coupling. Diag- 
onalization of the Hamilton matrix Ti, yields the sixfold 
spin-split valence band structure in the vicinity of the 
r point, depending on the magnetization orientation m 
and the strain in the (Ga,Mn)As layer. 

The m- and Ey-dependent normalized free-energy den- 
sity of the carrier system FM,c{eij , 'm) is obtained by first 
summing over all energy eigenvalues within the four spin- 
split heavy-hole and light-hole Fermi surfaces and then 
dividing the resulting energy density by M. The two 
split-off valence bands do not contribute to Fm,c because 
they lie energetically below the Fermi energy for common 
carrier concentrations. For biaxially strained (Ga,Mn)As 
layers grown pseudomorphically on (OOl)-oriented sub- 
strates, the tetragonal distortion of the crystal lattice 
along [001] can be fully described by the e^z component 
of the strain tensor using continuum mechanics. 

In order to compare the microscopic theory with the 
phenomenological description of the MA in Sec. IIII Al 
we consider the dependence of Fm on m with respect to 
the reference direction m.i-cf=[100]. Accordingly, we write 
the anisotropic part AFM.int of the intrinsic contribution 
-FM,int as 

An/,int = n/,c(e.., m) - Fm,c(£.., m = [100]). (7) 

In terms of the anisotropy parameters from Eq. ([2]), 
AF/v/jnt reads as 

Ai^M,int = B2^ml + B^\{ml + - I) + B^^mi. (8) 

For m rotated in the (001) and the (010) plane, Eq. ^ 
can be rewritten as 



AFM.inti^) = B^n (cos^ ip + siu'' ip-1) 



(9) 



and 



AFm i„t(6') = B2± cos^ e + S411 (sin'* 6-1) + cos'* 6, 

respectively, where we have introduced the azimuth 
angle if and the polar angle 9 with mx=sin9 cosip, 
TOj,=sin0sin(/?, and mz=cos9. We proceed by calculat- 
ing Ai^M,int numerically in the microscopic model as a 
function of (p and 9 with Ezz varied in the range —0.4% < 
Ezz < 0.3%, using typical values for p, Bq, and M. Equa- 
tions ^ and pU]) are then fitted to the resulting angular 
dependences using B2±, 54||i and B^^ as fit parameters. 
For the hole density we use the value p—3.5 x 10^° cm~'^, 
for the exchange-splitting parameter Bq=—2S meV, and 
for the magnetization fj,oM=4:0 mT. Inserted into Eq. 
the latter two values yield ApNof3=—1.8 eV, in good 
agreement with the parameters used in Ref. [^. 

The results of the microscopic calculations are depicted 
by the solid symbols in Figs, [^a) and[l](b). Since the 
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FIG. 1; AFAf.int calculated as a function of the magnetization 
orientation and the strain within the microscopic model (solid 
symbols) for M (a) in the (001) plane and (b) in the (010) 
plane, f and Q denote the azimuth and polar angles of M , 
respectively. The solid lines are least-squares fit curves using 
(a) Eq. dSl) and (b) Eq. JTOl) with B2X, B^\\, and B4X as fit 
parameters, (c) The anisotropy parameter _B2x obtained from 
the fit shows a pronounced linear dependence on Szz. 



variation of Ai^M,int('/') with Szz is found to be marginal, 
only one representative curve calculated with £2^^—0.4% 
is shown in Fig. [Ija). It clearly reflects the fourfold 
symmetry of Ai^jv/.int within the (001) plane. By con- 
trast, AFM,int{9) in Fig. [T](b) strongly depends on Szz- 
The solid lines in Figs. \V[a.) and \V[h) are least-squares 
fits to the calculated data using Eqs. ^ and pH)) . re- 
spectively. The perfect agreement between the micro- 
scopic results and the fit curves demonstrates that the 
intrinsic part of the free energy calculated within the 
microscopic theory can be well parameterized by B2±, 
B4U , and B4±. It thus justifies the phenomenological ap- 
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proach in Eq. taking into account only terms up to 
the fourth order in m. Whereas the cubic anisotropy pa- 
rameters i34||«i?4^«— 30 mT obtained from the fit are 
not substantially affected by Szz , the uniaxial parameter 
i32_L exhibits a pronounced linear dependence on Ezz, as 
shown in Fig. [ijc). The slope of B2^{£zz) is strongly 
influenced by the exchange-splitting parameter Bq and 
the hole density p, as will be discussed in more detail in 
Sec. UVCl Fig.Hb). 

The microscopic calculations show that for Szz < 
(tensile strain) Ai^M.int exhibits two equivalent minima 
for m oriented along [001] and [001], which become more 
pronounced with increasing tensile strain. In contrast, in 
the regime of compressive strain (szz > 0), the minima 
occur for m along [100], [100], [010], and [010]. Thus, 
the theoretical model correctly describes the well known 
experimental fact that for sufficiently high hole densities 
and low temperatures the magnetically hard axis along 
[001] in compressively strained layers turns into an easy 
axis in tensily strained layers. Note, however, that for a 
quantitative comparison between experiment and theory 
the extrinsic contributions to AFm have also to be taken 
into account. 



IV. RESULTS AND DISCUSSION 

In the following, the experimental data obtained for 
the (Ga,Mn)As samples under study are discussed. 



A. Lattice parameters and strain 

In order to study the structural properties of the 
(Ga,Mn)As/(In,Ga)As/GaAs samples, reciprocal space 
maps (RSM) of the asymmetric (224), (224), (224), and 
(224) reflections were recorded using HRXRD. Figure [2] 
exemplarily shows an RSM contour plot of the (224) re- 
flection for a nearly unstrained (Ga,Mn)As layer with 
£22=— 0.04%, depicting separate peaks for the GaAs sub- 
strate, the (In,Ga)As buffer, and the (Ga,Mn)As layer. 
From the peak positions {h,l) and their shifts {Ah,Al) 
relative to that of the substrate the lateral and vertical 
lattice parameters of (Ga,Mn)As and (In,Ga)As can be 
determined using the relations 



ail = as(l — Ah/h), aj_ ~ as(l — Al/l), 



(11) 



where denotes the lattice constant of the GaAs sub- 
strate, h and / are the coordinates in fc-space referring 
to the reciprocal lattice vectors of GaAs along the [100] 
and [001] directions, respectively. In Fig. [21 the peaks 
of (Ga,Mn)As and (In,Ga)As are centered at the same 
value of h, confirming that the (Ga,Mn)As layer has been 
grown lattice matched to the (In,Ga)As buffer. The lat- 
eral lattice parameters ay are therefore the same in the 
(In,Ga)As and (Ga,Mn)As layers. The shift Ah from the 
substrate peak at h=2 to lower values is due to the strain 
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FIG. 2: Reciprocal space map around the asymmetric 
(224) reflections of a nearly unstrained (Ga,Mn)As layer 
(e^z^— 0.04%) grown on relaxed (In,Ga)As/GaAs template. 
h and I are the coordinates in fc-space in units of the recipro- 
cal lattice vectors along [100] and [001] in the GaAs substrate, 
respectively. Ah and Al denote the shift of the peak position 
relative to that of the substrate. In the plot, and Al are 
only depicted for the (Ga,Mn)As layer. 



relaxation in the buffer layer (a|| >as). The lattice param- 
eters a|| of the (In,Ga)As templates in the as-grown and 
annealed samples are plotted against the In content in 
Fig.[3ja). Apparently, post- growth annealing had no sig- 
nificant influence on a\\, which linearly increases with the 
In content. 

For both (Ga,Mn)As and (In,Ga)As, the HRXRD mea- 
surements yielded different values of {h,l) for the (224) 
and (224) reflections, revealing a tilt of the lattice to- 
wards the [110] direction,^ The tilt of the (Ga,Mn)As 
layer originates from an equal tilt in the (In,Ga)As buffer 
pointing to an anisotropic relaxation of the (In,Ga)As 
templates, typically found in layers grown on vicinal 
substrateSi^ The samples under investigation, however, 
were grown on non-miscut (001) wafers. As can be seen in 
Fig. mb), the measured tilt angles tend to higher values 
with increasing In fraction. It should be emphasized that 
it is imperative to take the tilt into account when deter- 
mining the lattice parameters in order to avoid erroneous 
results. This can be done by inserting into Eq. (Ilip the 
averaged values of the peak positions and shifts obtained 
from the (224) and (224) reflections. For the angle- 
dependent magnetotransport measurements, the influ- 
ence of the tilt is negligible since the tilt angles observed 
are smaller than 0.06 degree. 

The relaxed lattice parameter a^ei of a biaxially 
strained layer on (OOl)-oriented substrate is obtained 
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and the (In,Ga)As layersi^ In Fig. HJ^a), Orei is shown 
for the (Ga,Mn)As layers as a function of the In con- 
tent in the (In,Ga)As buffer. It is found to be nearly 
unaffected by the (In,Ga)As template underneath. The 
fluctuations in the values of Crci are mainly attributed 
to slight variations of the growth temperature. For all 
(In,Ga)As templates under study, the degree of relax- 
ation defined by i?=(a|| — as)/(arci — Os) was above 80%. 
Once the vertical and relaxed lattice parameters of the 
(Ga,Mn)As layers are known, the vertical strain e^^ can 
be calculated from the relation 



(a_L — ftrcl)/ ftrcl- 



(13) 



In Fig. mb), Ezz is plotted against the In content. The 
slight decrease of Orei and e^z upon annealing is supposed 
to arise from the outdiffusion and/or rearrangement of 
the highly mobile Mni.iO->iiJ^iHiii 



FIG. 3: (a) Lateral lattice parameter a\\ and (b) tilt angle 
towards the [110] direction of the (In,Ga)As buffer layer, plot- 
ted against the In content. The solid line (as grown) and the 
dashed line (annealed) are regression lines. 
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FIG. 4: (a) Relaxed lattice parameter flroi and (b) vertical 
strain e^z of the (Ga,Mn) As layer plotted against the In con- 
tent in the (In,Ga)As buffer. Post-growth annealing only 
leads to a slight decrease of the values. 



from the relation 



Orel 



Cii + 2Ci2 



'II 



C\\ -I- 2Ci2 



(12) 



where Cn and Cyi are elastic stiffness constants. As 
an approximation, we use the values Cii=11.90x 10^° Pa 
and (712=5.34x101" Pa of GaAs for both the (Ga,Mn)As 



B. Hole density and Curie temperature 

Determination of the hole concentrations p in 
(Ga,Mn)As is complicated by a dominant anomalous con- 
tribution to the Hall effect proportional to the normal 
component of the magnetization M. (anomalous Hall 
effect). To overcome this problem, magnetotransport 
measurements were performed at high magnetic fields 
up to 14.5 T. Assuming the magnetization to be satu- 
rated perpendicular to the layer plane at magnetic fields 
IJ,oH>4 T, the measured transverse resistivity was fitted 
using the equation 



Ptr 



(H) = RofioH + cip\ong{H) 



C2P?onf 



(H) (14) 



for the ordinary and anomalous Hall effect with Rq, ci, 
and C2 as fit parameters. Here Ro=l/ep is the ordinary 
Hall coefficient and piong the measured field-dependent 
longitudinal resistivity. The second term on the right 
hand side arises from skew scattering"^^'"^^ and the third 
term from side jump scattering^ and/or Berry phase 
effectsi^ As mentioned in Sec. [TTl the Curie tempera- 
tures Tq were inferred from the peak positions of the 
temperature-dependent resistivities piong' ^^'^^ Consider- 
ing that the Tc values thus obtained generally differ from 
those determined by temperature-dependent magnetiza- 
tion measurements)^ we estimate an error margin of up 
to 20%. Similar to a^ei, neither p nor Tq are signifi- 
cantly influenced by the strain as shown in Figure^ The 
values scatter around pag=3.5xl0^° cm~'^ and Tc=Q5 K 
for the as-grown samples and Pann=5.8xl0^'' cm~'^ and 
Tc—91 K for the annealed samples. It is well known that 
the lattice constant, the hole density, and the Curie tem- 
perature strongly depend on the concentration of MuQa 
acceptors, Muj double donors, and other compensating 
defects such as AsQa antisites. The insensitivity of Oroi, 
p, and Tc to strain in the as-grown and annealed samples 
under study suggests the assumption that strain has no 
significant infiuence on the incorporation of Mnca, Mnj, 
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FIG. 5: Hole density p and Curie temperature Tc of 
(Ga,Mn)As plotted against the strain e^^ for the as-grown 
(solid symbols) and the annealed samples (open symbols). 
The fluctuations in the values are mainly attributed to slight 
variations of the growth temperature. For p we estimate an 
error margin of about ±10% and for Tq an error margin of 
up to ±20%. 



and AsGa- At least the sum of the changes caused by the 
different constituents seems to be unahered. Moreover, 
the insensitivity of Tq with respect to supports theo- 
retical predictions that the magnetic coupling should be 
unaffected by strain, since the corresponding deforma- 
tion energies are expected to be too small to significantly 
enhance or reduce the p-d kinetic exchange interaction.^"^ 



C. Anisotropy parameters 

Experimental values for the anisotropy parameters 
i3ooi=-B2±+-Bd, i?4||, i34j_, and i?iio were determined 
by means of angle-dependent magnetotransport measure- 
ments. A detailed description of the corresponding pro- 
cedure is given in the Refs. 21 and 12^ It can be briefly 
summarized as follows. The longitudinal and transverse 
resistivities piong and ptrans, respectively, are measured 
as a function of the magnetic field orientation at fixed 
field strengths of /io-ff=0.11, 0.26, and 0.65 T. At each 
field strength, H is rotated within three different crys- 
tallographic planes perpendicular to the directions n, j, 
and t, respectively. The corresponding configurations, 
labeled I, II, and III, are shown in Fig. [HI The vectors 




FIG. 6: The angular dependence of the resistivities was 
probed by rotating an external magnetic field H within the 
three different planes (a) perpendicular to n, (b) perpendic- 
ular to J, and (c) perpendicular to t. The corresponding 
configurations are referred to as I, II, and III. 



form a right-handed coordinate system, where j defines 
the current direction, n the surface normal, and t the 
transverse direction. For current directions j \\ [100] and 
j II [110], the resistivities can be written as^^ 

Piong = Po + Pirri^j + P2"^?l + Pa^^l + Pim^^ + pziu^jinl, 

(15) 

Ptians = P6m„ + pTnijUit + psml + pgrrijmtml, (16) 

where rrij, rrit, and m„ denote the components of m 
along j, t, and n, respectively. At sufficiently high mag- 
netic fields, the Zeeman energy in GM{m) dominates 
and the magnetization direction m follows the orien- 
tation h of the external field. The resistivity param- 
eters Pi (z=l,...,9) are then obtained from a fit of the 
Eqs. (fT5|) and (fT6|) to the experimental data recorded at 
0.65 T. With decreasing field strength, the infiuence of 
the MA increases and m more and more deviates from 
h. Controlled by the magnetic anisotropy parameters, 
the shape of the measured resistivity curves changes and 
Bqqi, i34j|, i?4j_, and -Buo are obtained from a fit to the 
data recorded at iJ,oH=0.26 and 0.11 T. In the fit pro- 
cedure, m is calculated for every given magnetic field 
H by numerically minimizing Gm with respect to m. 
Figure [7] exemplarily shows the angular dependence of 
Piong and ptrans for a nearly unstrained (Ga,Mn)As layer 
(e22=-0.04%) with H rotated in the (001) plane (con- 
figuration I) and j II [100]. The experimental data are 
depicted by red solid circles and the fits by black solid 
lines. 

Applying the procedure described above to the whole 
set of (Ga,Mn)As layers under study, the resistivity pa- 
rameters Pi (z=l,...,9) and the anisotropy parameters 
Booi=i?2±+-Bdj -S41IJ B4±, and Sjio were determined 
as functions of the vertical strain Ezz- The results for 
the resistivity parameters were extensively discussed in 
Ref. i22i . In the present work, we exclusively focus on the 
anisotropy parameters. 

Figure [5Ja) shows the values of the parameter 
BQ()i=B2±+Bd describing the uniaxial out-of-plane 
anisotropy. For both the as-grown and the annealed sam- 
ples, a pronounced linear dependence on Szz is found in 
qualitative agreement with the microscopic model calcu- 
lations presented in Fig. [TJc). For zero strain, the cubic 
symmetry requires -B2J.=0 and the extrinsic parameter 
i3d~60 mT is inferred from the intersections between the 
regression lines (dotted lines) and the vertical axis. If 
shape anisotropy was the only extrinsic contribution to 
the term of Fm, as assumed in Sec. IIII Al the value 
Bd~60 mT would correspond to a sample magnetization 
of /ioM==2i3d~120 mT. This value, however, exceeds the 
saturation magnetization determined by SQUID mea- 
surements by a factor of ~3. Figure [5] shows as an ex- 
ample the SQUID curves obtained for the as-grown and 
annealed samples with £22^0.2%. At the moment, the 
reason for this discrepancy is not yet understood. We sus- 
pect, however, that it might be related to the vertical gra- 
dient in the MA mentioned at the end of Sec. |TT1 Assum- 
ing the value of B^ to be nearly the same for all samples 
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FIG. 7: (Color online) Resistivities piong and ptrans recorded 
from a nearly unstrained (Ga,Mn)As layer witli Szz —— 0.04% 
at 4.2 K and j \\ [100] (red solid circles). The measurements 
were performed at fixed field strengths of iJ.oH—0.11, 0.26 
and 0.65 T with H rotated in the (001) plane corresponding 
to configuration I. The blacli solid lines are fits to the exper- 
imental data using Eqs. (I15p and (|16|l . and one single set of 
resistivity and anisotropy parameters. 



under consideration, the strain-dependent intrinsic pa- 
rameter B2± is obtained by subtracting from the mea- 
sured i^ooi data. In Fig. [8l[b), the values of i?2_L derived 
in this way are shown together with model calculations 
performed within the microscopic theory (see Sec lIIIBj) . 
Using for p the mean values Pag=3.5xl0^'' cm^'^ and 
Pj^jj„=5.8x10^° cm~^ (see Sec. lIVBt . the calculated val- 
ues of i?2± are found to be in quantitative agreement 
with the experimental data if Bq—~2S meV is chosen 
for the as-grown samples and i?G=— 39 meV for the an- 
nealed samples. The corresponding curves are depicted 
by the solid and dashed lines in Fig. ^h) . Comparing the 
exchange-splitting parameters with the respective hole 
densities, we find 



^Cann 



Pa 



Pa 



1.7. 



(17) 



For the annealed samples, theoretical results are also 
shown for —32 mcV and —23 meV, demonstrating that 
the slope of B2±{szz) is drastically reduced with decreas- 
ing Bq. Remarkably, normalization of the experimen- 
tally derived values of B2± to the corresponding hole con- 
centration p yields the same linear dependence of B2±/p 
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FIG. 8: (a) Dependence of the uniaxial out-of-plane 
anisotropy parameter i?ooi=S2±+Bd on Szz for the as- 
grown (sohd circles) and annealed samples (open circles). 
i3d=60 mT is inferred from the intersections between the cor- 
responding regression lines (dotted lines) and the vertical 
axis, (b) Anisotropy parameter B2± obtained by subtract- 
ing 60 mT from the measured Booi data. The lines are model 
calculations for B2± performed within the microscopic theory 
described in Sec. IIII 61 using values for the parameter Bq as 
shown and the averaged p values from Fig. [5] (c) Anisotropy 
parameter B2_l=Booi — Bd normalized to the hole density p. 
Experimentally, the same linear dependence of B2±/p on Ezz 
is obtained for both the as-grown and the annealed samples. 
The dotted Une represents a linear regression. 



on Ezz for both the as-grown and the annealed samples 
as shown in Fig.[5]Jc). We thus find, at least for the range 
of hole densities and strain under consideration, the ex- 
perimental relationship 



B2± = Kpsz 



(18) 



with K—\.b7x\Q~'^'^ Tcm"^. If we assume a linear rela- 
tion BQ—ApfiM/dgfiB oc p, in accordance with Eq. (fTT)) . 
Eq. (jl8p can be reproduced by the microscopic theory. 
Note however, that this is not trivial, since Ai^Af.int ex- 
plicitly depends on both Bq and p. Anyhow, the relation 
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FIG. 9: SQUID curves of the as-grown and annealed samples 
with ezz«0.2%, measured at 5 K for H oriented along [001]. 
The saturation magnetization of ~40 mT is representative for 
the whole set of (Ga,Mn)As samples. The values of the co- 
ercive fields, derived from the hysteresis loops (not shown), 
were found to be below 10 mT. 
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FIG. 10: Dependence of the cubic anisotropy parameters (a) 
S4II and (b) B4± on Ezz for the as-grown (solid circles) and 
annealed samples (open circles). The lines are model calcula- 
tions within the microscopic theory using the same values for 
p and Bq as in Fig. [8jb). 



Bq oc p demands a future detailed investigation. 

The experimental values of the fourth-order parame- 
ters B^ and B4± are presented in Figs. [TUTa) and fTPf b). 
respectively. Whereas i?4|| only slightly varies between 
—40 and — lOmT, B4± exhibits positive values close to 
10 mT for e^z^— 0.15% and negative values of about 
-20 mT for e^^>-0.15%. The lines depicted in Fig. [TO] 
were calculated using the same values for p and Bq as 
in Fig. [SJ^b) . Obviously, the change of sign of B4± does 



FIG. 11: Dependence of the uniaxial in-plane parameter -Biio 
for the as-grown (solid circles) and annealed samples (open 
circles) . 



not appear in the theoretical curves. Apart from this dis- 
agreement, the experimental data of the as-grown sam- 
ples are again well reproduced for Bq=—23 meV. In the 
case of the annealed samples, now the splitting param- 
eter Bg=—32 meV yields a much better description of 
the measured data than Bq——39 meV, found to fit the 
strain dependence of i?2±- 

In view of the perfect quantitative interpretation of 
B2± by the mean-field Zener model, this small discrep- 
ancy and the fact that the change of sign of B4± is not 
reproduced by the calculations, should not be overesti- 
mated. The experimentally observed change of sign of 
B4± may be caused by extrinsic influences, not accessible 
by the model, such as the increasing density of threading 
dislocations in the (Ga,Mn)As/(ln,Ga)As layers with in- 
creasing In concentration. Moreover, we cannot rule out 
that the change of sign is an artifact of the experimental 
method for determining the relatively small anisotropy 
parameter B4±. The same problem has also been re- 
ported in Ref. 19, where the authors did not extract 
information on the fourth-order out-of-plane anisotropy 
parameters from their ferromagnetic resonance measure- 
ments, since the corresponding contributions to the MA 
were masked by the much larger contributions of the 
second-order out-of-plane term. As to the differing val- 
ues of —39 meV and —32 meV for Bq in the case of the 
annealed samples, it should be pointed out that for high 
hole densities the Fermi energy shifts deep into the va- 
lence band. Therefore, the values of the free energy ob- 
tained by using a 6 x 6 fc • p effective Hamiltonian become 
increasingly unreliable when analyzing higher order con- 
tributions to the free energy. One should not jump to the 
conclusion that this has to be interpreted as a deficiency 
of the Zener model itself. 

In agreement with the data presented in Refs. [l^, l4ll . 
and I42I the experimental values of the extrinsic uni- 
axial in-plane parameter Bjiq are much smaller than 
those of the cubic in-plane parameter B411 at 4.2 K. 
As shown in Fig. 1111 they scatter around zero with 
-10mT<Bno<10mT. 
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Due to the strong dependence of -Bqoi on Szz-, the out- 
of-plane axis [001] becomes magnetically harder with e^z 
increasing from —0.4% to 0.2%. Neglecting the small 
influence of Siio in our samples at 4.2 K, the critical 
strain e™*, where a reorientation of the easy axis from 
out-of-plane to in-plane occurs, can be estimated from 
the condition _Booi+S4j_=i34|| . We obtain for the as- 
grown (Ga,Mn)As layers under study e™*=— 0.13% and 
for the annealed layers e™*=— 0.07%. Remarkably, these 
values are very close to the e^z values in Fig. [10] where 
i34_L changes sign. 

Since the MA sensitively depends on the individual 
growth conditions, care has to be taken when compar- 
ing the values of anisotropy parameters published by 
different groups. Keeping this restriction in mind, the 
data presented in this work are in reasonable agreement, 
e.g., with the results obtained by Liu et alJ^ for a rep- 
resentative pair of compressively and tensily strained 
(Ga,Mn)As samples with 3% Mn. 



D. Magnetostriction constant 

Changing the magnetization of a ferromagnet, e.g. by 
an external magnetic field, leads to a variation of its ge- 
ometrical shape. For crystals with cubic symmetry, the 
relative elongation A in a given direction /3 can be ex- 
pressed in terms of the magnetization orientation m and 
the magnetostriction constants Aioo and Am along [100] 
and [111], respectively, according toA^^ 



+3A111 {m^iriyP^Py + niymzPyPz + m^mzPxPi;}^) 



Starting from Eq. (|T8)) . we are able to determine Aioo 
defined hy^ 



Ainn — — - 



9 C12 — Cii 



(20) 



where ai denotes the magnetoelastic coupling con- 
stant. First-order expansion of the free-energy density 
F{eij,m) with respect to Eij and continuum mechanics 
yield the relation B2A_M=ezzai{l + C11/2C12). Thus, 
Eq. (pni) can be rewritten as 



Ai 



2KMp 



9(Ci2-Cii)(l + Cii/2Ci2) 



(21) 



Inserting the experimental values for p, /ioAf, and if, 
we obtain Aioo~^3 ppm for the as-grown samples and 
Aioo~— 5 ppm for the annealed samples. In Ref. O . 
we already deduced an approximately constant value of 
Aiii«5 ppm below 40 K decreasing to zero at higher 



temperatures (40 K<T<Tc«85 K) via applying piezo 
stress along the [110] direction of a piezoelectric ac- 
tuator/(Ga,Mn)As hybrid structure. Our results for 
both magnetostriction constants are close to the values 
Aioo=— 11-3 ppm and Aiii=8.1 ppm reported by Mas- 
manidis et al.— 

V. SUMMARY 

A series of (Ga,Mn) As layers with 5% Mn was grown on 
relaxed graded (In,Ga)As/GaAs templates with In con- 
tents up to 12%. In this way, the vertical strain Szz in the 
(Ga,Mn)As layers could be gradually varied over a wide 
range from —0.38% (tensile strain) to 0.22% (compres- 
sive strain). The strain was found to have no significant 
influence on the hole concentration, the Curie tempera- 
ture, and the relaxed lattice parameter. Angle-dependent 
magnetotransport measurements were performed to de- 
termine the uniaxial and cubic anisotropy parameters. 
B2± turned out to be proportional to both the strain and 
the hole concentration. From this linear dependence, the 
magnetostriction constant Aioo was determined. While 
-B211 and -B411 are nearly strain independent, Bi± changes 
sign when the magnetic easy axis flips from in-plane to 
out-of-plane. 

Microscopic calculations of the free-energy density 
were performed based on the mean-field Zener model of 
Dietl et al.^ They justify the approximations made in 
the parameterization of the free energy by considering 
only terms up to the fourth order. The strain-dependent 
anisotropy parameters derived from the calculations were 
found to be in good quantitative agreement with the ex- 
perimental results. In the case of the as-grown samples, 
the values of B2± and i34|| are even perfectly reproduced. 
The quantitative comparison between the comprehensive 
set of experimental data and the microscopic calculations 
may be considered a valuable contribution to the ongo- 
ing controversy on impurity band versus valence band in 
(Ga,Mn)As4^iil 

Using the orientation of the magnetization as the ba- 
sic information bit of a non-volatile memory, the tay- 
loring and manipulation of the MA is of special impor- 
tance. As shown in this work, the choice of an appropri- 
ate (In,Ga)As template allows for an adjustment of Ezz 
close to the critical value e™*. Then, the magnetization 
direction can be switched from in-plane to out-of-plane or 
vice versa, e.g., by minute variation of the piezo stress in 
a piezoelectric actuator/ (Ga,Mn) As hybrid structure^. 
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